A simple analytical description for the cross-tie domain wall structure. 
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A closed form analytical expression for the magnetization vector distribution within the cross-tie 
domain wall in an isotropic ferromagnetic thin film is given. The expression minimizes the exchange 
energy functional exactly, and the magnetostatic energy by means of an adjustable parameter (wall 
width). The equilibrium value of the wall width and the film thickness corresponding to the tran- 
sition between the Neel and the cross-tie walls are calculated. The results are compared with the 
recent experiments and are in good qualitative agreement. 



The structure of magnetic domain walls (transition regions separating magnetic domains) was in a focus of extensive 
research in 1930s - 60s and during the time a lot of useful results on the static structures of various wall types as 
well as their dynamic properties were obtained. Many of the results on the theory of the domain walls can be found 
in books @, In particular, there are analytical expressions for the distribution of the magnetization vector inside 
of the one-dimensional Bloch and Neel domain walls, which are the starting point for calculating their static and 
dynamic properties. 

The cross-tie domain wall was also observed in a number of experiments a that time (see Ref. [| and references 
' therein) and also with modern high resolution techniques [5], Jj), and are usual in thin and ultra-thin magnetic films 
important for modern applications However, there is no j|| (also see p. 163 in Ref. |2|) closed form analytical 
expression for the cross-tie wall structure. This is, probably, due to the fact that magnetization distribution even in 
straight cross-tie domain wall is two-dimensional (unlike one-dimensional ones of Bloch and Neel walls). It means, 
the structure of such a wall is defined by a system of non-linear integral (due to the long-range dipolar interactions) 
partial differential equations, and there is no way to reduce this system to a single equation (as it was done for the 
Neel wall |§). 

Consider a thin film having the thickness h made of soft (isotropic) magnetic material, in the Cartesian coordinate 
system X, Y, Z chosen in such a way that 0Z axis is perpendicular to the film plane. The parameters of the material 
entering the calculation are the exchange constant C, and the saturation magnetization constant Mg. If the film is 
thin enough, so that h is of the order of a few exchange lengths Le = \J C/Mg, the dependence of the magnetization 
distribution on Z can be neglected and the task becomes essentially two-dimensional. If we also neglect the dipolar 
interaction for a while, the magnetization distribution M(r) = Mgin(r), \m\ = 1, r = {X, Y} is defined by the 
. minimum of the exchange energy functional: 

i=x,y,z 

where the integration runs over all X — Y plane. The last expression is given using the parametrization of the 
magnetization vector field by a complex function w(z, z) of a complex variable z = X+iY, i = \f— T, line over a variable 
means complex conjugation, so that m x +im y — 2w / (1+ww) and m z = (1— ww) / (1+ww) , d/dz = (8/dX—id/dX)/2, 
d/dz = (d/dX + id/dY)/2. It is easy to see that the energy functional ([!]) is scale- invariant, and the exchange energy 
of the magnetization distribution rh{r/c) or w(z/c,~z/c), where c is real constant, is independent on c. If we include 
the dipolar interaction back into the picture this scale-invariance breaks. The dipolar interaction will be considered 
. in this paper to determine the scale of the magnetization distribution obtained by minimizing (0) . 
The system of Euler equations for the function w(z, z) giving extremum to the functional (H) is 



d ( dw\ 2w dw dw 



M ', dz \ dz ) 1 + ww dz dz 

It is easy to see that any analytic (in a sense that the Cauchy-Riemann conditions dw/dz — are satisfied) 
function of a complex variable satisfies the equation (||). The partial solution of this class was first written as a 
rational polynomial in z by Belavin and Polyakov [^| and can be thought as a superposition of vortices. There are 
other solutions whose applications to the magnetism of small cylindrical particles are given elsewhere [fl0|. 
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FIG. 1: Distribution of the magnetization vector within the cross-tie domain wall having with c = 1. 



Consider the function 

w c-t( z ) — i tan(z/c), (3) 

while the Cauchy-Riemann conditions are satisfied for this function, and, therefore, it minimizes the func- 
tional (0) exactly, it is not a Belavin-Polyakov (BP) soliton because it can not be represented as a rational 
polynomial in z and, consequently, has infinite energy when integrated over the whole space (while the en- 
ergy of BP solitons is always finite). Even though the total exchange energy of the magnetization distribution 
wc-t is infinite, it is extremal, because wq-t satisfies (0). Also, the corresponding fa — {m x ,m y ,m z } 
{-tanh(2y/c),sin(2X/c)/cosh(2y/c),cos(2X/c)/cosh(2y/cj| (compare ]ll)) plotted in Fig. | has the topology of 
the magnetization distribution in the cross-tie domain wall 
The exchange energy per unit of the wall area is finite 



7 CX = ^— = I dx I dY * = M 2 hX 2^z ( 4 ) 
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where dimensionless variables A = Ls/h and £ = cit/h were introduced. This energy term (|4|) forces the wall to 
expand infinitely as it has the smallest value only for c, Q — > oo. 

The evaluation of the magnetostatic energy is straightforward, but a little bit voluminous, it was done by calculating 
the density of magnetic charges (having both volume and surface contributions) , solving the Poisson equation for the 
scalar magnetic potential in the Fourier representation and performing the convolution of the density of charges and 
the potential to obtain the energy. The final result can be represented as 

y* = M i h Ko, (5> 

7T J Jo z; 2 cosh 2 (7ru/2)(l + u 2 + v 2 ) 

One of integrals in /(C) can be taken analytically, the other is rapidly converging and easy to take numerically. The 
representation (^) was chosen for its compactness. It is also possible to build a simple analytical approximation 
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Bloch domain wall 




FIG. 2: The normalized energy per unit wall area (solid lines) of the equilibrium cross-tie, Neel and Bloch domain walls as 
a function of film thickness. The dashed line shows the normalized width of the cross-tie wall. The dotted lines show the 
equilibrium cross-tie domain wall energy and width calculated using the approximate magnetostatic function. 



for /(C) 



/(C) = CxC + C 2 (l-e- 2 ^)C 2 
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(1 + u 2 )(l + cosh7ru) 

1 



o 4(l + u 2 ) 3 / 2 cosri (ttu/2) 



0.17753... 
0.12119. 



(7) 
(8) 

(9) 



which is asymptotically correct for C <C 1 and produces the worst case error at C — ► oo around 5%. 

The total energy per unit wall area 7 = j cx + r y dlp was minimized numerically with respect to the parameter C, and 
the result is shown in Fig. || in comparison to the energy per unit area of the Neel and Bloch domain walls calculated 
according to the model of Dietze and Thomas |l2|] . As it could be expected for films thicker than a certain value 
h w 2.3Le the cross-tie domain wall takes over the Neel one as the lowest energy configuration. Another interesting 
fact is that in the shown region of film thicknesses (where roughly the assumption of uniformity of the magnetization 
distribution along the thickness, dM/dZ = 0, can be expected to hold) the energy of the cross-tie domain wall is 
lower than the energy of Bloch wall. This suggests that the transition from the cross-tie wall to the Bloch one takes 
place at larger thicknesses, where dM/dZ = can not be safely assumed anymore. 

As it is visible from Fig. ^ the approximate expression for magnetostatic function (f?j) after the energy minimization 
produces results almost indistinguishable from the exact ones and having only a slight deviations for extremely small 
film thicknesses. Thus, the expression (^) is safe to use instead of ([|) for most practical estimations of the cross-tie 
wall magnetostatic energy. 

The qualitative validity of the presented model is supported by the good fit of the image signal (roughly proportional 
to the magnetization vector component parallel to the wall plane) presented in Fig. 10 of Ref. [| to the empirically 
chosen taxih.(Y/6 w ) function. In fact, m x (X,Y) = — tanh(2Y"/c) exactly for the solution (^). Thus, the cross-tie wall 
width, as defined in Ref. 0, is S w = c/2. The permalloy film in that experiment had roughly h/ Le = 60nm/18.2nm 
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3. Reading from Fig. || we get the theoretical width of the cross-tie wall S w = 4.82v£:/(27r) w 14nm, which is of the 
order of the measured value j| of 42 ± lOnm. It is expected that better quantitative correspondence can be obtained 
for thinner films and also by including the anisotropy energy (which is small for permalloy) into the minimization. 

The presented simple model for the magnetization distribution in the cross-tie domain wall can be a starting point 
for analytical estimations of its static and dynamic properties, especially for short wavelength excitations on the 
cross-tie wall background, where the fact that (||) minimizes the exchange energy functional exactly can lead to 
significant simplifications. It can be also useful for interpretation of MFM images, and for supplying initial conditions 
for numerical finite-element micromagnetic simulations. 
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